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Abstract 

We consider a generalization of the Heath Jarrow Morton model for the term struc- 
ture of interest rates where the forward rate is driven by Paretian fluctuations. We 
derive a generalization of Ito's lemma for the calculation of a differential of a Pare- 
tian stochastic variable and use it to derive a Stochastic Differential Equation for the 
discounted bond price. We show that it is not possible to choose the parameters of 
the model to ensure absence of drift of the discounted bond price. Then we consider 
a Continuous Time Random Walk with jumps driven by Paretian random variables 
and we derive the large time scaling limit of the jump probability distribution func- 
tion (pdf). We show that under certain conditions defined in text the large time 
scaling limit of the jump pdf in the Fourier domain is uit(k, t) ~ exp {— &/ (ln(fct)) 2 } 
and is different from the case of a random walk with Gaussian fluctuations. 
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1 Introduction 



The Arbitrage Pricing Theory (APT) is a generalization of the Capital Asset 
Pricing Model (CAPM), and it states that capital assets are priced according 
to a set of factors. It does not say what these factors are, just that only 
idiosyncratic (firm specific) risk is priced. So, two different stocks could have 
two different APT formulations. 

However, the law of one price states that two identical assets should have the 

same value on the financial market and deviations from rational behav ior of 

th I 

traders for example, significantly larger returns on Fridays the 13 . Lucev 



(2000) rave been evidenced. Yet, market efficiency is an assumed prerequisite 



for C APM to hold and cont e mpor ary financ e still uses APT as a quant itative 
tool. Heath- Jarrow-Mortonl (^392) (see also Musiela- Rutkowskil (|l99lh have 



formulated a general theory (to be referred to as the HJM model) and a 
unifying framework for the existing APTs. 

Presuming the forward interest rate and the bond price to be a continuous 
time stochastic process described by a Stochastic Differential Equation (SDE) 
they derive the generic form of the SDE that ensures that the bond price, 
related to the zero coupon bond (risk-less security issued by the central bank 
at a fixed interest rate called the spot rate), has a mean zero as a function of 
time. 



The HJM model has been generalized by Sornettei ( 1998f ). who, drawing on 



some ideas from string theory, derives a general condition for a class of second 
order SDEs that must be satisfied to ensure absence of arbitrage opportunities. 
He obtains a whole class of, nonlinear second order SDEs that, in some cases, 
can be solved exactly and that als o provide extensions to the theory of pricing 



options M usiela- Rutkowsk i (1991) 



These models have a rich structure but are based on Gaussian fluctuations. In 
this paper we develop a corresponding theory that allows the fluctuations to 
have an arbitrary distribution, in particular they may have power-law behavior 
at the high end. 



2 Theoretical analysis 



The purpose of this section is to provide overall understanding of stochas- 
tic processes that involve asymptotically large, power-law distributed, fluc- 
tuations. The analysis lacks mathematical rigor and draws more on physical 
ideas of Langevin processes (even though they are not explicitly mentioned in 
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the text) than on the abstract theory of stochastic processes. The derivation 
of strict mathematical foundations for the analysis (e.g. foundations for the 
Lemma and the Statement in the next section) is left for future work. 



2. 1 Definitions and heuristic reasoning 



The forward interest rate f(t,T) is an interest rate that is agreed at time t 
and that gives the trader a right to trade in an asset at time T (to be termed 
as maturity). The underlying assumption in the HJM model is to describe 
f(t; T) through a following SDE: 

t+e t+e 

fit + e; T) = f[t; T) + J a (9;T) d9 + J a (9; T) V {6)d9 (1) 

t t 



where e > 0. Here the functions a and er, called drift and volatility are in 
general dependent on the process itself, and the fluctuations 5n =V i6)d6 have 
a zero mean, are uncorrelated with one another for different times and conform 
to a generic, non-Gaussian distribution with a density piSn). An asset which 
at maturity has a value one Euro has at some earlier time t a smaller value 
that is expressed through a discount factor involving the forward rate (1) as 
follows: 

B(t;T)=B(t,f(t,T);T)=ex V l- J /(*;£) j# (2) 



The bond price is related to the price of a zero-coupon bond M(t) whose 
stochastic dynamics is given by a SDE that involves the spot rate of interest 
fid, 9) and some fluctuations: 

t t t 
M(t) = exp \ - f f(9, 9)d9 + f V (9)d9 \ (3) 



and conditions are derived for the discounted bond price M(t)5 (t; T) to have 
a zero drift in order to ensure absence of arbitrage opportunities. The zero- 
drift condition is also referred to as the "martingale property" and states that 
the conditional mean of the bond price at time s > t discounted back to time 
zero equals the discounted bond price at time t. Here E [M(s)5(s;T) |JF t ] = 
M(t)i?(t;T) and the conditioning is on the so called "filtration" Tt of infor- 
mation meaning a s et of information about the bond price up to time t (see 
iMusiela-Rutkowskil (|l99lh . 
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The main problem is now to derive a SDE for the bond price (2). For comput- 
in g a differen tial of a function of Gaussian fluctuations one can use the lemma 
bvlitc1(|l997h. This lemma states that the differential is obtained by expanding 
the function into a Taylor series to the first order in time and to the second 
order in the fluctuations and by replacing differentials of the fluctuations via 
their averages. If the function of the stochastic process f has a form given 
by equ ation (2) we may use a theorem by Girsanov [see iMusiela-Rutkowski 
that simplifies the computation. 



If, however, the fluctuations are not Gaussian, in particular if the likelihood 
for them to be large is much larger than in the Gaussian case then the higher 
moments and auto-covariances of the process may be infinite (as in the case of a 
Levy or a Pareto distribution) and the procedure for computing the differential 
may fail. It may also happen that the probability distribution of the squared 
fluctuation is of the same kind as that of the fluctuation itself which will 
certainly change the structure of the SDE that describes the bond price (2). 
In the following we examine the mathematical difficulties that arise. Since 
we are not aware of any generalization of Ito's lemma that can be used in a 
straightforward manner for the calculation of the differential of the bond price, 
we must perform the computations from scratch. The procedure consists in 
first working out a small change of the bond price, secondly in expressing it 
via the forward rate and finally in performing the limit e — > . Here the change 
reads : 



B{t + t- T) — B(t;T) = B(t;T) 



e r +e /(^K exp [_y" 5e/(t . e)de j _! 



(4) 



where 5 t f(t; £) := f(t + e; £) — f(t; £). Now we expand the right hand side in 
a Taylor series and retain terms proportional to e neglecting higher powers of 
e. The expression in square brackets on the right hand side of (4) reads: 



- l-iy ' 1 



J f(t + 6, o« + E -^f- - / <W; 0# 

t+e T t+e T t+e 

J f(t + e,Od£- J J a(9;0ffldS- J J a (9; £) V (9)d9d£ 
t t t t t 

+ E^r f// <T{e,t)v((f)d0dt\ (5) 



vt t 



Now the whole problem consists in a proper treatment of the products of 
fluctuations Yli=i i n the last term on the right hand side in (5). This 

is the point where the analysis differs essentially from the Gaussian case. We 
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assume the fluctuations to conform to a Levy (or Paretian) distribution. We 
set 5r] = Levy^(0, 1, 0)dt which means that V— 5r]/dt is (parameters in brack- 
ets from the left to the right) a Paretian distribution with a characteristic 
exponent /x, skewness zero, dispersion one and mean zero (to be termed stan- 
dard Levy variable). The Paretian distribution has a pdf whose asymptotic 
expansion reads: 

1 00 7r dt n 
p( 6 v) = - £sin(-/m)r(/m+ 1)-^ (6) 



From this expansion follows an important property for the product of many 
Paretian variables. 

Lemma: 

(A) If Sr, = Levy^O, 1, 0)dt then (5 V ) n = Levy^jO, 1, 0)(dt/n) 

(B) For independent Levy variables 5r]i and 5i] 2 we have 

Z = ^ 2 ~Levy M (0,l,0)(^) 2 . 

Both statements regard the high end of the distribution. 
Proof: 

(A) From the conservation of elementary probabilities /y (?/)<% = fx{x)dx it 
follows that the pdf of the n-th power Y = X n of a Paretian variable 
X = Levy M (0, 1, 0)dt reads: 

i ) = fx ^ = f*(y 1/n ) = 1 y a ^ dtm = 1 y a ^ dfm 

^y^^Levy^ /n (0,l,0) (7) 

where a m denote coefficients in the expansion (6) for the pdf of a Paretian 
variable. 

(B) If Xi,X 2 = Levy M (0, l,0)dt are independent Paretian variables then the 
pdf of their product Z = X X X 2 reads: 



z z 

fz(z) = J J S(z - x 1 x 2 )fx(x 1 )fx(x 2 )dx 1 dx2 
1 1 



z 



fdx! z r dx 1 (aidt) (a-^dt) {a-^dty 

= J — f ^ )fX ^ = J — ^ (g/flM+l = 1^ HZ) ^ 

^Z~(rft) 2 Levy M (0,l,0) (9) 

where (8) holds only for large values of z since we have replaced the 
pdfs by the first term of the asymptotic expansion (6). From the last 
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expression in (8) it follows that ,up to a logarithmic correction, Z is a 
Paretian with a squared dispersion q.e.d. 

The Lemma implies that the only factors proportional to the first power of e 
in the last sum on the right hand side of equation (5) are those which contain 
"diagonal" terms of the kind (5f](6i)) n . Indeed from (B) it follows that the non- 
diagonal terms, ie those including products (Si](6i)) ni (Si](62)) n ' 2 . . . {Si](6 p )) np 
of differentials of r\ evaluated at different values of 6 will all be proportional to 
higher powers of e (since they are products of independent Paretian random 
variables) and can, therefore, be neglected. We emphasize the difference to the 
Gaussian case: 



Statement The n-th power of a Gaussian with variance dt is a random vari- 
able with variance dt n . The nth power of a Paretian with a dispersion dt also 
has a dispersion dt. 



We now obtain, after taking the limit e — > the following SDE for the bond 
price: 



d t B(t,T) r/ , } . ^ „ ~ (-1)™ 1 1 



71=1 



/(*, t)-j a{t; i)di + exp - lj a{t, £R % it) 



- 1 (10) 



In the next step we carry out the same procedure for the discount factor M(t) 
defined in (3). Expanding M(t + e) — M(t) in a Taylor series to arbitrary 
order in the variable /„* (f{6,0)- V {6)dd) d9, as in (5), we notice that after 
expanding the n-th (with n > 2) power of the integral we get some terms that 
contain f(9,9) and other terms that contain only the fluctuations V (6). The 
former terms will be of order e 2 and may be neglected. The only term left then 
is the nth power of the integral over fluctuations. This we rewrite as a sum 
of increments 8r)(t + idt) at times (t + idt). This corresponds to a division of 
nodes of the interval [t,t + e] into e/dt bins. Thus we have: 



( / v (9)der = n 

t i =1 



e/dt 

Y,Sr}{t + ijdt) 



e/dt 



J2(Ht + Jdt)) n + 0(e 2 ) = (V it)fe + 0(e 2 )(ll) 



i=i 



Now notice that on the grounds of Lemma and/or of the Statement that 
the power terms (5rj(t+jdt)) n are dt/nLevy^/ n (0, 1, 0) whereas all other terms 
5rj(t + j\dt) ni • . . . • 5rj(t + j p dt) np that involve p > 2 different time values are, 
since they are products of p independent Paretian variables, on the grounds 
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of Lemma (B) of order of dt p . Since their number is a polynomial of order 
(e/dt) p the total will be of order e p and can, therefore, be neglected. 

In the limit e — > only the powers of Paretian variables survive and we arrive 
at the following SDE for the discount factor: 



d t M(t) 
M{t) 



-/(M)+exp(^ (t)) 



:i2) 



Adding the logarithmic derivatives (12) and (10) we obtain the logarithmic 
derivative of the discounted bond price V(t) = M(t)B(t;T), ie the price at 
time zero of a bond that at maturity T pays off one Euro and has, at some 
intermediate time t, the value B(t; T). 



d t V{t) 
V(t) 



-2 + exp 



}<r(t,£)dA V, (t) 



exp(% (*))(13) 



This is an SDE for the discounted bond price that we must solve and then 
require the Martingale condition E [V(s) — V(t) \Ft\ = for every s > t to 
be satisfied in order to ensure non-arbitrage opportunities in the market with 
non Gaussian fluctuations. Note that the time evolution of the logarithm of 
V(t) is driven by two competing factors namely a negative "deterministic" 
drift corresponding to the first two terms on the right hand side in (13) and a 
positive fluctuations (the last two terms). Therefore the Martingale condition 
boils down to requiring the average fluctuation to be equal to the deterministic 
drift. Note also that, due to an exponential enhancement of the underlying 
fluctuations (£), the fluctuations in (13) are very large; and huge when 
compared to Gaussian noise. It is, therefore, not clear whether the absence of 
arbitrage opportunities can be indeed ensured (see item (B) in the discussion 
section 2.2). 



2.2 Connections to other works in this subject matter 



At this stage we make some remarks regarding the connection (A),(B) of our 
result to works by other authors and to an abstract representation (C) of 
infinitely divisible distributions. 

(A) The equatio ns (10) and (13) cannot be directly related to the result (12) 
page 154 in iKleinerd ((20041 ) . This is because the author assumes the 
logarithmic characteristic function <p(k) := logE exp(ik V (t)) of the 
fluctuations T) (t) is an analytic function in k or, in other words, that all 
moments of the fluctuations are finite ((f) (t)) p ) < oo, (see expansion of 
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the "Hamiltonian" - equation (6) page 153 in iKleinertl (J200J)) whereas 
we have <p(k) = logE zfcLevy^O, 1,0) = — which is not analytic as 
k — ► (some of the constants ci, c 2 , . . . in expansion (6) page 153 may be 
infinite) . 

(B) iBieleckil f)2004h . using the framework of the abstract probability theory, 
considers a Gaussian HJM model with several forward interest rates, so 
called default-able gi(t;T) i = 1,...,K rates and a default free rate 
f(t; T), subject to an ordering condition gx{t\ t) > . . . > g±(t; t) > f(t; t) 
on the respective spot rates of interest. They define a discounted Default- 
able Corporate Bond (DCB) Z(t; T) that equals the consecutive dis- 
counted default-able bond prices Zj(i;T). This is related to the default- 
able forward rates of interest gi(t; T), and the discounted default free rate 
of interest /(i; T) in time intervals 7$ <t < r i+ i\ intervals that are deter- 
mined by the condition that the differences between the default-able bond 
prices Di(t\T) and the default free B(t;T) price do not exceed a given 
threshold. In other words they consider a process of "credit migrations" 
where the DCB is taken as D^(t; T) at time inception and is subsequently 
replaced by Di(t; T) with i < K once differences benchmarked on the de- 
fault free price B(t; T) exceed a given threshold (see Fig.l ). In addition, 
certain rates 5i G [0, 1] of recovery (returns of losses from the Treasury 
Bond - central bank) are assumed in each time interval t G [Tj, r i+1 ]. Fi- 
nally conditions on the absence of arbitrage opportunities for the DCB 
are deduced. 



Default Free Bond 
1st Defaultable Bond 
2nd Defaultable Bond 
Corporate Bond 




-20 ■ 



Fig. 1. Time evolution of the logarithm of price of the Discounted Corporate Bond as 
considered in lBieleck] ( 20041 ) . The evolution starts with the most highly default-able 
bond (Diamonds) and goes down to the less default-able (Triangles) and default free 
(Boxes) bonds once a threshold difference from the default free logarithm of price 
is reached. 
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This is certainly a very sophisticated "real world" model yet it does not 
compare to our approach because in the presence of large Paretian fluctu- 
ations (see Statement), that are presumed in our approach, absence of 
arbitrage opportunities cannot be ensured. Indeed, a closer explanation 
of the time evolution (13) of the discounted bond price V(t) ensures that 
it is impossible to choose the drift a(t; T) and the volatility a(t\ T) such 
that the discounted bond price as a whole has a zero drift, ie the mean 
of a small increment V(s) — V(t) for s > t of the discounted bond price 
conditioned on the past history of the price process equals zero under a 
certain probability measure. This follows from a following "waving hand 
argument" . Due to an exponential enhancement of, already very strong, 
Paretian fluctuations; the mean of the right hand side of equation (13) 
is infinite. Therefore the positive "stochastic jumps" dictated by the last 
two terms in equation (13) are, for sufficiently large times t, arbitrarily 
large and cannot be offseted by any finite value of the negative "deter- 
ministic" drift coefficient — J t T a(t; £)d£ — 2. Thus the discounted bond 
price drifts toward positive values like a biased random walk. Since it is 
impossible to damp the biased movement of the discounted bond price 
V(t). Further theoretical analysis of the model will be pursued in a dif- 
ferent direction namely we shall work out a distribution of random times 
needed for the biased movement to reach a certain threshold j: value of 
the log-price, known as first hitting times 



T y = min [t \ln(V(t) = p] 



(14) 



of the price. This can be done by writing down a Fokker-Planck d ifferen- 
tial equation, an equation that may include fractional derivatives Samkol 
( 19931 ): Meerschaert ( 2001 ) and that describes the spatio-temporal evolu- 
tion of the probability for the discounted bond price to have a given value 
at a given time (to be termed occupation probability). The equation may 
be solved in the asymptotic limit of large times an d the pr o babili ty dis- 
tribution of the first hitting times will be derived Rednerl (j2001l )] from 
it. This approach will be pursued in section (3). 
(C) Since the stochastic process (1) that describes the forward rate of interest 
is an Infinitely Divisible Distribution (IID) it must have a representation 
in terms of the Ley y-Khintchin formula (e.g. Theorem 3.1.11 page 41 in 
Meerschaert ( 200ll) ) for the log-characteristic function of an IID. Indeed 
the stochastic dynamics (1) of the forward rate can be considered as a 
a compound Poisson process f(t;T) = a(t;T) + cr(t;T)J2jli^Vj with 
Nt = Poisson(c e ) and 5r]j = Levy„(0, 1,0) with a very large rate c e . 
This corresponds to taking a small positive e > and a measure (f)(dx) 
in the Levy-Khintchin formula equal to <f)(dx) ~ l/x a+1 . In this case 
c e — dx/x a+1 can be arbitrarily large for e sufficiently small. Thus the 
jump pdf p(Sr)) ~ c~ l l/ (5r]) a+1 is related to a Paretian distribution for 
large values of the jump size Srj. 
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3 Occupational probability of the discounted bond price 



Now we formulate an equation for the occupational probability V (x, t) = 
P (In V(£) = x) for the logarithm In V(£) of the discounted bond price in (13) 
to attain a value x at time t. Note that since the process lnV(i) is a Con- 
tinuous Time Random Walk (CTRW) with jumps occurring at random times 
the probability density function (pdf) <fi n (x,t) of reaching x at time t after 
having performed n steps is a convolution of the pdf (p n -i (x',f) and a jump 
probability uj t (x — x',t — t'). Here 



oo t 



<Pn (X,t) 



uj t (x — x' , t — t') 4> n -i (a/, t') dt'dx' 



(15) 



oo 



Now we use the fact that the occupational probability V (x, t) is a convolution 
in time of the probability density of reaching x at time if after an arbitrary 
number of steps n and th e probabi l ity de n sity of b e ing un moved in the time 
interval [*',*]■ We obtain |Montrolll (|l965h : iBalescul (|l997h ] the following for- 
mula for the Fourier-Laplace transform P(k,s) := T x L t [V (x,t)} (k,s) (with 
respect to x and to time t respectively) of the occupational probability. 



P(k,s) 



1 - 



uj t (k, s) 



(16) 



where u t (k,s) = J r x C t [u t (x,t)} = ff^ J °° exp(ikx — st)u t (x,t)dxdt is the 
Fourier-Laplace transform of the jump probability density and <ft(s) = c^(0, s) 
is the Laplace transform of the waiting time pdf. Thus in order to find the 
occupational probability we need to work out the Fourier-Laplace transform 
of the jump pdf. 

Denoting E := J 4 T a(t,^)d^ and A := a(t; £)d£ we express the jump proba- 
bility as: 



uj t (x -x',t- t') = P (In V(t) - In V(t') = x - x') 




exp 



S V, (Ol + exp(^ (0)) d£ = x-x'+ (A(0 + 2K 



P £ («P [-SO] + exp(0)) =(x-x' + j (A(£) + 2)d^j j 

°r °r ( M N ( r ' 

/•••/ * f E(^' + e^)- + /(A(0+2K 



A? 



(17) 
(18) 

(19) 



n p(Cs)m 

3=1 
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where the equality (17) = (18) follows from the SDE (13). The next equality fol- 
lows from the division of the time interval [t', t] into N intervals of length At/N 
and from the definition of a Riemann integral. The final equality (19) = (20), in 
which 5 means the Dirac delta function, follows from the fact that the values 
of fluctuations Q =rj^ (£,•) at intermediate times — t' + j(t — t')/N are iid 
standard Levy variables. Performing the Fourier transform of the jump pdf 
with respect to x — x' results in replacing the delta function in the integrand 
in (20) by an exponential exp {ik . . .} and produces a following path integral 



oo oo ( t 

X] (21) 



oo oo f t \ 

Q t (k,t- = — J . . . J exp \ikj (OK D 



— oo — oo V V 



where the "Hamiltonian" H(^,V^ (0) reads: 



= (exp [-E V, (0] + exp(^ M (£))) - (A(£) + 2) 



(22) 
(23) 



Here we divided the "Hamiltonian" into two parts the "stochastic" part 
Tis^lfx (£)) that depends on the fluctuations ?7 M (£) and the "deterministic" 

part Hd(£) that depends on the drift function A(£)- The measure D 



is 



expressed via pdfs of standard Levy variables (£,•) as follows: D 



lim 



Hj=iP(Vti. (€j))d ^1 n (0) Note tliat since tlie "Hamiltonian" depends 
only on the fluctuations (0 and not on their time derivative (0 1 the 
path integral (21) is expressed via ordinary Riemann integrals only. Here: 



u> t (M-0 



2vr 



N^oo 



f 



N 



(24) 



where 

oo oo 

f(fc):= | exp{^ s (C)}p(CR = / exp{^(e^ + e- s O}p(C)rfC (25) 

— oo — oo 

and ro(iV) is a number of sub intervals into which we split the interval [t' , t] so 
that the the path integral does not diverge when N — > oo. 

1 This would correspond to a propagator of the Schrodinger equation in quantum 
mechanics with a Hamiltonian that only depends on the coordinates of the particle 
and not on its velocities. 
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In order to get additional insight into the solution of the problem in our further 
calculations we replace the integrand in the definition of f(k) by a sum of two 
terms in which one of the terms is neglected. The applicability of this 
approximation is yet to be examined. 

Approximation: 



f(fc)~-(S(i,*) + 3(-E,fc)) 



(26) 



where 



oo oo 

d(a,k) := J exp {tke<} p(()d( = ^- J e xpH£n(-)^r(-</aK(27) 



In the second equality on the right hand side in (27) we have substituted for 

and used the integral representation 
p(0 = 1/(2tt) ff^,exp(-z£ - IfHdf of the Levy distribution. We then have 
changed integration limits in one of the integrals by making use of the Cauchy 
theorem from complex analysis and replaced one of the integrals by the Gamma 
function making use of its integral representation. The minus sign in the last 
expression in (27) stems from changing the integration limits after applying 
the Cauchy theorem. In further calculations we exploit the following results: 



Laplace Transforms: 

A [t a ] (s) = T(a + l)/s a+1 



Ct 



'ikt 



(s) = l/(s-ik) 



C t [f(kt)\ (s) = 1/kCt [f(t)\ (s/k) 
1 

S (ln(s))%t^o V ^ (Hs)) p 



(ln(f))» 



is) 



' i \n oo p 



-at 



(\n(kt)Y 



1 



1 



E 



c 



(«+«) re) B siN)) p 



(28) 
(29) 
(30) 

(31) 
(32) 



The Gamma Function 



r(x)r(l - x) = 7r/sin(7rx) for xeC 



(33) 



Convolutions of Laplace transforms 
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(i«i)=i (34) 
\s sj s 

( s " kmRf) = ' (35) 

(5(111(5))™ s(ln(s)) m J = s(ln(s))"+ m J Dm ' P (ln(s))P ' (36) 

where m p := E£=i(-lW +1 , n P := E~i andD m , p :=(m + p- 1)1/ (r. 

Here <g> denotes a convolution 

(/(s) ® g( s ))( s ) := Jq f(£)g(s — £)d£ and the remaining parameters are defined 
and the results are derived in Appendix. 

Now we work out the Laplace transform of the function f(t). We get: 

f(s) = Ct W)\ (s) = \ (5(1, s) + £(-E, S )) (37) 



where 



1 00 1 

~ — 1 f ITT 

5 <°- S > = 2^ / ^Hfl" WMflp^™^ (38) 



—OO 



^ 7exp(Hy + V(2 n )0 Sin( l / ;' I ' / ( ? ^ (39) 
2as J sinh(7r£/a) 



-00 
00 



± J exp (-|f r + «/(2a)0 - — i— ■ ^fy^ di (40) 
2as J smh(7r£/a) 4 



-00 



In the first equality (38) we make use of (28) and of the result (33) concerning a 
product of Gamma functions. In (39) we write the term s^/ a as exp(?^/aln(s)) 
and we disregard the imaginary part owing to the fact that the result is real. 
In (40) we separate the integrand into a product of terms that will be handled 
later on. 

Now we find a small- s approximation to the function $(a, s). Since for small 
positive s the last part of the integrand in (40) is a very "high frequency" 
periodic function of £ with a maximum (ln(s)/a) at £ = and a period 
of length 27ra/|m(s)| we apply the following approximations (to be termed 
steepest descent approximations) for our calculations. 

b(l) Expand the first term := exp (— + ir/ (2a)£) in the integrand (40) in 
a Taylor series in £ to the second order around zero. Note that the maximum 
of the whole integrand (the product of all three terms in (40)) is bigger than 
zero but, since the expression sin(£(ln(s)/a))/£ behaves for small positive s 
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as 6(£/a), the maximum of the whole integrand will be arbitrarily close to 
zero for sufficiently small s. 

b(2) Confine the integration range to the first period [— ira/\ ln(s)|, na/\ ln(s)|] of 
the "high-frequency" periodic function in the last term of the integrand. 

b(3) Replace the middle and the last terms in the integrand by a small £ expan- 
sion around £ = 0. 



This yields a following result: 



7ra/| ln(s)| 

las J \ l J 7r a 

— 7ra/| ln(s)| 



1 ln(s) 



2as 7r 



7ra 



ln(s) 



1 7T 2 . . . . 9 1 1 . / 1 

- - —| 3 (0+)|a 2 - ; Na + -O 



6 



s ln(s) 5 



ln(s) z 



(41) 



This calculation is incomplete in that we did not explain sufficiently the na- 
ture of approximation used by computation of the integral. In fact we should 
have expanded the integrand around its maximum £ = £o( s )> calculated the 
integral and only then considered the result in the limit of small positive s. 
This would only change the value of the parameter J (0+) in the result. In 
our calculations this parameter is left unknown. 

Now we invert the Laplace transform f(s) and, on the grounds of (31) and of 
(37), we get: 

ii'f^-l^^'-Sf^fc) (42) 



where & is some parameter that is unknown. Inserting (42) into (24) we obtain 
the jump probability in the limit of large times: 



Co t (k,t) = ^-e- lk ^ nd{ ^ -Mm 



2tt 



2tt J (In 



N 



kt -jy 



■w(N) 



1 e - lfe />(0+2K . exp 



8. 



2tt 



(43) 

(44) 
(45) 



14 



where in (45) we have chosen the number of subintervals zu(N) as w(N) = [kt)^ +l 
to ensure that the limit is finite. Now we use result (32) to calculate the 
Fourier-Laplace transform Co t {k, s) of the jump pdf and to obtain the Fourier- 
Laplace transform of the occupational probability P(k, s) of the discounted 
bond price from (16). We assume for simplicity that the function A(£) does 
not depend on time. 



u t (k, s) = Ct 



e lkVt ■ exp 



n=0 



-ikVt 



(\n(kt)) 2n 
1 



(ikV + s) 



1 -8. 



(\n(kt)) 2 



(ln(!)) 2 



+ 




(46) 
(47) 
(48) 



where V : = (A + 2). Inserting k = in (48) we get <j)(s) = a> t (0, s) — 1/s what 
implies that the waiting time pdf (f>(t) is constant <f>(t) = 1 as a function of 
time and is not normalisable / °° <p(t)dt = oo! Because of this constructing the 
occupational probability V (x, t) from (16) in a process where the transition of 
the walker to x in time t consists of performing a certain number of jumps and 
of being immobile for some time may be impossible. This is because the Fourier 
Laplace transform u t (k, s) may be bigger than one and thus the formula (16) 
cannot be applied. 

Note that if the drift A(£) satisfies a condition: A + 2 = then the Fourier 
transform of the jump pdf decays as u t (k, t) ~ exp {— &/ (\n(kt)) 2 } This decay 
is much slower than the inverse of any positive power of t. The interpretation 
of this result and fitting the model to high-frequency financial data is left for 
future work. 



4 Master equation for the discounted bond price 



In this section we derive a partial differential equation that is satisfied by the 
jump pdf ujt (x, t). The derivation bases on the representation (21) of the jump 
pdf as a path integral. Note that equation (24) can be written as follows: 



&t(k,t - t') = — exp 



2tt 



ik / (A(o+2)de 



• exp 



h {m) \k)di 



(49) 
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where 

^\k) :=-ln(f(Jfe)) = In 



exp {tk(e ( + e" EC )} p(()d( 



(50) 



The function ^{k) can be viewed as a "log-characteristic function" related 
to a certain integral transformation of the jump pdf p(C), a transformation 
defined by equation (25). If we now assume A(£) and £(£) to be constant as 
functions of £ then from (49) we obtain an equation for the Fourier transform 
of the jump pdf 

d t u t {k, t) = (A + 2) (—ik)u> t (k, t) - fjV\k)u> t {k, t) (51) 

and, after taking the inverse Fourier transform, an equation for the jump pdf 
itself: 

dtut{x,t) = -{A + 2)d x uH{x,t)-S)P\-id x )u t {x,t) (52) 



where the operator id x ) is defined by expanding the function ^j^(k) in 

a Taylor series around k = and by replacing k by —id T . The ma ster equation 
(52) corresponds to formula (112) on page 17 in Kleinertl (|2002[ ). 



In future work we will check if the master equation (52) includes fractional 
derivatives. In particular we will in vestigate if it ca n be related to di s tributed 
order fr actional diffu s ion eq u ations C lechkinl |2002); Gorenflol ( 20031 ); Sokolov 
( 2004 ): iMeerschaertl (|2004h : ISchefflerl ifeool that are widely studied due to 
their connection to ultraslow diffusion. 



5 Conclusions 



We derived the large time limit of a Continuous Time Random Walk where 
the motion of the walker is driven by a stochastic process with Paretian fluc- 
tuations. We applied the CTRW model to the Heath- Jarrow-Morton model 
for the term structure of interest rates and showed that if the forward rate 
of interest is driven by Paretian fluctuations then the bond price follows a 
stochastic process that behaves like a biased random walk. This means that in 
our model arbitrage opportunities cannot be canceled out by an appropriate 
choice of the drift and the volatility parameters of the HJM model. Further 
work will be devoted to deriving asymptotics of distributions of first hitting 
times in the CTRW model and to developing a theory of pricing options in 
the Black and Scholes framework. 
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At the end of writing the manuscript we became aware of a different work [ 
Eberlein Raiblel l|l999h ] devoted to the subject matter. We will refer to that 



paper in our future work. 
6 Appendix 



In this section we derive mathematical results regarding the small-s expansion 
(36) for a convolution of Laplace transforms, and Laplace transforms (31), (32) 
of a power of an inverse logarithm and of a function that involves exponentials 
and powers of inverse logarithms. 

At first we compute the convolution that is used for computing the Laplace 
transform of the path integral (21). Here: 



ln(s) ln(s) 

= / = 1 / rfi£(n)( e *--ta(»h (54) 

J z m (s - e z ){ln(s - e z )) n s{\n(s)) n J z m 



oo 



:-d- t ^,(£1) (55) 



s(ln(s)) n 



where in (53) = (54) we substituted for ln(£), then we introduced an auxiliary 
function £,^(x) : = (1 - x)-\l + ln(l - x)/\n(s))' n and finaly in (54) = (55) 
we substited for —z. Now we expand the function £( n ) in a Taylor series 
Z^ n \x) = X^Lo lp(n)x p , insert the expansion into (55) and we get: 



l) m ^P m -\(n) 7 _ m 

i w mexp^-wjdw (56) 



s(M*)) n £i sP 

r — pln(s) 

We recall that the integrals (56) are integral representations of the truncated 
Gamma function 

oo 

T{x,m) := / w~ m exp(— w)dw (57) 



which has a following asymptotic expansion: 
p/ \ , ^vv ^. (m + g-l)! 1 

r(x,m) = exp(-x)}J-l) g (m _ 1) , xm+(? (5* 
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Inserting the expansion (58) into (56) we obtain the result (36) q.e.d. Note 
that the coefficients l p (n) in the expansion of the function also depend 
on s. They have a form: l p (n) = ^ p n ^ (l/log(s)) where are polynomials 
of order p and ty p n ^ = 1. 

Now we derive the Laplace transform (32) of a function that involves an ex- 
ponential and powers of the inverse natural logarithm. We calculate: 



3 — at 



(\n(kt)) n 



I \ 1 ( *r- I a + S \ o-+ I a + S \ 



(59) 



where 



J n (a, a) 



exp(— at) 



and X+(a,a) := I ^ ~f rft (60) 
nV 7 y (ln(t)) n 



In the following we compute the integrals I+\a,a) and Z^(a:, a). For com- 
puting the integral T+\a,a) we consider a following auxiliary integral: 

/exD f — ott) 
(\n(t)) n tP dt such that J n (« ? a ) = ^n,o(a,a) (61) 



The integral J n , P (oi, a) differentiated with respect to the parameter p yields 
the truncated Gamma function (57). We have 



d n J n ,p{a, a) 
dp n 



I 



,-at (-Ht)T dt 

(ln(t)) n tP 



■l) n a p ~ l Y{aa,p) = (-l) r 



a 



(62) 

Erl^r^ 1 "'"' (63) 



q=0 



(ao)« (p-1)! 



where on the right hand side in (62) we used the identity l/t p = exp(— p\n(t)) 
and in (63) we used the asymptotic expansion (58) of the truncated Gamma 
function. Integrating (63) n times with respect to p from p to infinity and 
using the fact (61) that d m J nyP (a, a)/dp m \ p=oo = for m = 0, . . . , n we get: 



J n>p (a,a) = (-l) r 



a 



E 



■i) 4 



o ( aa ) 



pip) 



q p 



^ in(a) n(e+«),e 

s=0 



(64) 
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where 



oo oo oo 



p 6 p 



We insert the last expression on the right hand side in (65) into (64), we 
evaluate the integral over £ by successive integration by parts and we obtain: 



rr f \ (-1)" e— ^ (-!)« Wj\ q ( P ) 



where 



Wr^) ■■= e- l li(Z + s) (67) 

s=0 



and W^li qip) denotes the /th derivative evaluted at £ = p. 

Now we set p — 0. Since the polynomial W n _i )g (£) has (67) a leading power 

£ n_1 then the derivatives W^ l l l q (0) for I — 0, . . . , n — 2 at zero are equal to 

zero. This allows to transform the expression (66) for p = to the following 

form: 



J n ,o(a,a) = Z+(a,a) = 
1 



-V 



VV n+q-l l U J ^ -gg 



n-1)! a (Ma)) n S (Ma))' (««) ? 



(68) 



It is easy to work out the form of the coefficients W^ l _^\o). They take the 
form 



in 



+ /-l)!(g-l)!^_" 1 1) for 1>1 
(n- l)!^ for 1 = 



(69) 



where 

£ for /> land #°> = 1 (70) 

l<£i<-<&<r 



Inserting (69) into (68) and changing the summation index q = I + r to r we 
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get: 



2+(a, a) = (-!)"- 



a (ln(a)) r 



E 



ra + Z - 1 1 A(aa) 



I 



(Ha)) 1 



(71) 



where 

A (a) := 



E 

r=0 



(72) 



Now we calculate the integral X n (a, a). The procedure is similar to that used 
in the first part of the appendix for deriving result (36). The integral reads: 



X n (a, a) = (-1)" J exp(-ae z )— dz = 

- In (a) 

("If E^T-(^ + ^""^(-(Z + 1) ln(a), n) = 

( j fro « 1 j (n-l)! (-(Z + l)ln(a))^ 

(ln(a))«^(/ + l)! ^ (n-l)! ((Z + 1) ln(a))« 1 J 

1 ~ ( w + f_i)! j ~ (_ a )g q^+i 
(R^F h (^-1)! (In(a))' ^ g! (g + 1)' +1 



(73) 

(74) 
= (75) 



(77) 



where in (73) we substituted for — ln(t), in (74) we expanded the exponential 
exp(— ae~ z ) in a series of e~ z and transformed the resulting integrals to the 
truncated Gamma function, in (75) we made use of the asymptotic expansion 
(58) of the truncated Gamma function and finaly in (76) and in (77) we 
collected together terms with the same power of the inverse logarithm. 



X n {a, a) 



a(m(a))' 



E 



n + l-l) Bi{aa) 



I 



(Ha)) 1 



(78) 



where 



i+i 



(79) 
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The proof of the last equality in (79) is left to the reader. Having derived the 
results (71) and (78) we complete the computation of the Laplace transform 
(59). The result obviously depends on the parameter a that divides the in- 
tegration domain into two parts. Since we are interested in the small s limit 
of the Laplace transform we choose a to be a decreasing function of s, e.g. 
a = k/s. This conforms to the intuition since a has to be large enough for the 
asymptotic series in (72) to converge. We obtain the following result: 



e 



1 1 1 x> (n + l-1 \ (Bi(0) + (-l) n e- e Ai{6)) 



(\n(kt)) n 



1 1 00 

(a + s) (ln(f))»§ 



(m(f)) p 



(80) 



where 9 = (a + s)/s. Therefore the proof of result (32) is finished. The Laplace 
transform (31) of a power of an inverse logarithm follows immediately by 
inserting k = 1 and a = into (80). The coefficients C p and 6 P in (32) and 
(31) are expressed via coefficients B p (8) and A p {9) as follows: 



C p = (n + p- l)\/(p\{n - 1)!) • (B p (9) + (-l) B e-%(0)) (81) 
t p = (n + p - l)!/(p!(n - 1)!) ■ (B p (l) + (-1) V 1 ^!)) (82) 
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